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Abstract 

On the basis of frequentist analyses of experimental constraints from electroweak pre- 
cision data, {g — 2)fi, B physics and cosmological data, we predict the masses of Higgs 
bosons and SUSY particles of the constrained MSSM (CMSSM) with universal soft 
supersymmetry-breaking mass parameters, and a model with common non-universal 
Higgs masses (NUHMl). In the CMSSM we find preferences for sparticle masses that 
are relatively light. In the NUHMl the best-fit values for many sparticle masses are even 
slightly smaller, but with greater uncertainties. We find that at the 95% C.L. all colored 
particles are in the reach of the LHC. While the light Higgs boson is bounded from 
above by ^125 GeV, the heavy Higgs bosons could well escape the LHC searches, 
but might be accessible at the ILC. 
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INTRODUCTION 



Supersymmetry (SUSY) [[ESIl] is one of the favored ideas for physics beyond the Stan- 
dard Model (SM) that may soon be explored at the Large Hadron Collider (LHC). In sev- 
eral recent papers f^J3, 6], we presented results from frequentist analyses of the param- 
eter spaces of the constrained minimal supersymmetric extension of the Standard Model 
(CMSSM) — in which the soft supersymmetry-breaking scalar and gaugino masses are 
each constrained to universal values mo and mi respectively (see [6] for a comprehen- 
sive list of references) — and the NUHMl — in which the soft supersymmetry-breaking 
contributions to the Higgs masses are allowed a different but common value (see [6] for 
a comprehensive list of references). Other statistical analyses in these models can be 
found in ^ S, §, [M JIMMMS^ and Markov Chain Monte Carlo 

(MCMC) analyses in flllSHMElSi^^ 

Here we review the results presented in They include the parameters of the best- 
fit points in the CMSSM and the NUHMl, as well as the 68 and 95% C.L. regions ap- 
plyingthe phenomenological, experimental and cosmological constraints. These include 
precision electroweak data, the anomalous magnetic moment of the muon, {g — 2)p_, B- 
physics observables (the rates for BR{b sy) and BR(5„ TV^), mixing, and the 
upper limit on BR(5i )), the bound on the lightest MSSM Higgs boson mass, 

Mfj, and the cold dark matter (CDM) density inferred from astrophysical and cosmolog- 
ical data, assuming that this is dominated by the relic density of the lightest neutralino, 
n^h^. In [5] we also discussed the sensitivities of the areas of the preferred regions to 
changes in the ways in which the major constraints are implemented. We found that 
the smallest sensitivity was to the CDM density, and the greatest sensitivity was that to 



DESCRIPTION OF OUR FREQUENTIST APPROACH 



We define a global 'y} likelihood function, which combines all theoretical predictions 
with experimental constraints: 

M ( fobs _ rfit n2 

+ ;c2(suSY search limits) + V ^^^^^^-^^ (1) 

/ Cy(/SM,)^ 

Here A'^ is the number of observables studied, C, represents an experimentally measured 
value (constraint) and each Pi defines a prediction for the corresponding constraint that 
depends on the supersymmetric parameters. The experimental uncertainty, a(C,), of 
each measurement is taken to be both statistically and systematically independent of 
the corresponding theoretical uncertainty, o{Pi), in its prediction. We denote by X^i^h) 
and ;(^^(BR(5i l^t^)) the contributions from the two measurements for which only 
one-sided bounds are available so far, as discussed below. Furthermore we include the 
lower limits from the direct searches for SUSY particles at LEP [.36] as one-sided limits, 
denoted by ";t:^(SUSY search limits)" in eq. ©• 

We stress that in [6] (as in |0,[5l]) the three SM parameters /sm = {Acthad^'^^^z} are 
included as fit parameters and allowed to vary with their current experimental resolutions 
^(/sm) • We do not include as a fit parameter, which would have only a minor impact 
on the analysis. 

Formulating the fit in this fashion has the advantage that the probability, 
P{X^,Niiof), properly accounts for the number of degrees of freedom, A^dof^ in the 
fit and thus represents a quantitative and meaningful measure for the "goodness-of-fit." 
In previous studies PiX^i^dof) has been verified to have a flat distribution, thus 
yielding a reliable estimate of the confidence level for any particular point in parameter 
space. Further, an important aspect of the formulation is that all model parameters are 
varied simultaneously in our MCMC sampling, and care is exercised to fully explore 
the multi-dimensional space, including possible interdependencies between parameters. 
All confidence levels for selected model parameters are performed by scanning over 
the desired parameters while minimizing the X^ function with respect to all other 
model parameters. The function values where X^i^) found to be equal to Xmin~^^X^ 
determine the confidence level contour. For two-dimensional parameter scans we use 
Ax^ = 2.28(5.99) to determine the 68%(95%) confidence level contours. Only ex- 
perimental constraints are imposed when deriving confidence level contours, without 
any arbitrary or direct constraints placed on model parameters themselves. This leads 
to robust and statistically meaningful estimates of the total 68% and 95% confidence 
levels, which may be composed of multiple separated contours. 

The experimental constraints used in our analyses are listed in Table 1 in ['6']. One 
important comment concerns our implementation of the LEP constraint on M^. The 
value quoted in the Table, Mh > 114.4 GeV, was derived within the SM [|37|l . and is 
applicable to the CMSSM, in which the relevant Higgs couplings are very similar to 



those in the SM [|38L 1391] . so that the SM exclusion results can be used, supplemented 



with an additional theoretical uncertainty: we evaluate the X'^i^h) contribution within 
the CMSSM using the formula 



(l.lGeV)2 + (1.5GeV)2 



(2) 



with M)™" = 1 15.0 GeV for Mh < 1 15.0 GeV. Larger masses do not receive a X^i^h) 
contribution. We use 115.0 GeV so as to incorporate a conservative consideration of 
experimental systematic effects. The 1.5 GeV in the denominator corresponds to a 
convolution of the likelihood function with a Gaussian function, <I>i 5(jc), normalized 
to unity and centered around M^, whose width is 1.5 GeV, representing the theory 
uncertainty on M/^ [40]. In this way, a theoretical uncertainty of up to 3 GeV is assigned 
for ~ 95% of all M/, values corresponding to one CMSSM parameter point. The 1 . 1 GeV 
term in the denominator corresponds to a parametrization of the CL^ curve given in the 
final SM LEP Higgs result [l37|]. 

Within the NUHMl the situation is somewhat more involved, since, for instance, a 
strong suppression of the ZZh coupling can occur, invalidating the SM exclusion bounds. 
In order to find a more reliable 95% C.L. exclusion limit for in the case that the SM 
limit cannot be applied, we use the following procedure. The main exclusion bound from 
LEP searches comes from the channel e+e ZH,H bb. The Higgs boson mass limit 
in this channel is given as a function of the ZZH coupling in [41]. A reduction in the ZZh 
coupling in the NUHMl relative to its SM value can be translated into a lower limit on 
the lightest NUHMl Higgs mass, M^™"'*^, shifted to lower values with respect to the SM 
limit of 1 14.4 GeV. (The actual number is obtained using the code HiggsBounds [|42n 
that incorporates the LEP (and Tevatron) limits on neutral Higgs boson searches.) For 
values of M/^ < 86 GeV the reduction of the ZZh couplings required to evade the LEP 
bounds becomes very strong, and we add a brick-wall contribution to the x^ function 
below this value (which has no influence on our results). Finally, eq. (O is used with 
^hmit _ ^imit,o _|_Q^ ensure a smooth transition to the SM case, see [6] for 

more details. 

The numerical evaluation of the frequentist likelihood function using the constraints 
has been performed with the MasterCode which includes the following 

theoretical codes. For the RGE running of the soft SUSY-breaking parameters, it uses 
Soft SUSY p3|] . which is combined consistently with the codes used for the various 
low-energy observables. At the electroweak scale we have included various codes: 
FeynHi ggs M(x l44l l45l . l46|] is used for the evaluation of the Higgs masses and a^^^^ 
(see also [47. 48Ll49l l5(^. For the difference between the SM and the experimental value 



we used = (30.2 ±8.8) x lO^^" [51] based on e+e 



data, see also [SSIHl. A 
[El- Using 



new evaluation, including new BaBar data, yields Aa^ = (24.6 ± 8.0) x 10 
this value could have a small impact on our results. We note that recently a new T based 
analysis has appeared fs^, which yields a ~ 1.9 a deviation from the SM prediction. For 
flavor-related observables we use SuFla ['sv', '58] as well as Super I so [|59L [60ll. and 
for the electroweak precision data we have included a code based on [|6il[6SFinally,for 



dark-matter-related observables, MicrOMEGAs [[631 |64L |65|] and DarkSUSYj66 



have been used. We made extensive use of the SUSY Les Houches Accord [[68 
the combination of the various codes within the MasterCode. 



m in 



RESULTS FOR SPARTICLE MASSES 



For the parameters of the best-fit CMSSM point we find niQ = 60 GeV, = 3 10 GeV, 
Ao = 130 GeV, tan/3 = 1 1 and /i = 400 GeV, yielding the overall X^/^dof = 20.6/19 
(36% probability) and nominally Mf, = 1 14.2 GeV. The corresponding parameters of the 
best-fit NUHMl point are mo = 150 GeV, mi/2 = 270 GeV, Aq = -1300 GeV, tan/3 = 

11 and m^^ = = — 1.2 x 10^ GeV^ or, equivalently, /i = 1140 GeV, yielding = 
18.4 (corresponding to a similar fit probability to the CMSSM) and M/, = 120.7 GeV. 

We now review the results for the predictions of sparticles masses in the CMSSM and 
the NUHMl, which are summarized in Fig. [TJ The results for the CMSSM spectrum 
are shown in the left plot, and for the NUHMl in the right plot. We start our discussion 
with the gluino mass, nig. In both the CMSSM and the NUHMl, the best-fit points have 
relatively low values of nig ~ 750 and ~ 600 GeV, respectively. These favored values are 
well within the range even of the early operations of the LHC with reduced centre-of- 
mass energy and limited luminosity. However, even quite large values of nig < 2.5 TeV 
are allowed at the 3-a (Ax^ = 9) level (not shown in Fig. [T])- The LHC should be 
able to discover a glui no with nig ~ 2.5 TeV with 100/fb of integrated luminosity at 
^/s = 14 TeV H eI, and the proposed SLHC luminosity upgrade to 1000/fb of 
integrated luminosity at ^/s = 14 TeV should permit the discovery of a gluino with 
mg ~ 3 TeV [72]. However, Fig.[T]does demonstrate that, whilst there are good prospects 
for discovering SUSY in early LHC running [5], this cannot be 'guaranteed'. 

The central values of the masses of the supersymmetric partners of the u,d,s,c,b 
quarks are slightly lighter than the gluino, as seen in Fig. [B The difference between 
the gluino and the squark masses is sensitive primarily to mq. The reason is that the 
preferred regions of the parameter space in both the CMSSM and the NUHMl are in 
the ^I'-slepton coannihilation region [|5|,|6|] where wq < Here mo makes only small 
contributions to the central values of the squark masses. The SUSY partners of the left- 
handed components of the four lightest quarks, the qi, are predicted to be slightly heavier 
than the corresponding right-handed squarks, qn, as seen by comparing the mass ranges 
in Fig. [U As in the case of the gluino, squark masses up to ~ 2.5 TeV are allowed at 
the 3-a level. Comparing the left and right panels, we see that the squarks are predicted 
to be somewhat lighter in the NUHMl than in the CMSSM, but this difference is small 
compared with the widths of the corresponding likelihood functions. 

Turning now to the likelihood functions for the mass of the lighter stop, m^j , we find 
that it is shifted to values somewhat lower than for the other squark flavors. It can also 
be seen that the 2-a range of its likelihood function differ from those of the gluino 
and the other squarks, reflecting the importance of scalar top mixing. We recall that this 
depends strongly on the trilinear soft SUSY-breaking parameter A/ and the Higgs mixing 
parameter /i, as well as on the precise value of m/. 

In the case of the lighter stau fi, see its range in Fig. [H the mass is very similar to 
that of the LSP Xi the coannihilation region, but this is not the case in the rapid- 
annihilation H,A funnel region, see [6] for details. In the case of the NUHMl rapid 
annihilation is possible also for low tan/3, leading to larger values of mo than in the 
CMSSM also for relatively small values of m^j . 

The scalar taus as well as the other scalar leptons are expected to be relatively light. 
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FIGURE 1. Spectra in the CMSSM (left) and the NUHMl (right) The vertical solid lines indicate 
the best-fit values, the horizontal solid lines are the 68% C.L. ranges, and the horizontal dashed lines are 
the 95% C.L. ranges for the indicated mass parameters. 



as can be seen in Fig. [T] They would partially be in the reach of the ILC(500) (i.e. with 
y/s = 500 GeV) and at the 95% C.L. nearly all be in the reach of the ILC(IOOO) ItbIItI. 
This also holds for the two lighter neutralinos and the light chargino. 



PREDICTION OF HIGGS BOSON MASSES 



In Fig. [2] we display the favored regions in the (MA,tan/3) planes for the CMSSM 
and NUHMl. We see that they are broadly similar, with little correlation between the 
two parameters. Concerning tan/3, one can observe that while the best fit values lie at 
tan/3 ^ 1 1, the 68 (95)% C.L. areas reach up 30 tanjS ^ 30(50-60). The existing Higgs 



discovery analyses (performed in the various benchmark scenarios [|13L |14| . |75L udA ) 



cannot directly be applied to the (MA,tan/3) planes in Fig. [21 In order to assess the 



prospects for discovering heavy Higgs bosons at the LHC in this context, we follow the 
analysis in [77], which assumed 30 or 60 fb^^ collected with the CMS detector. For 
evaluating the Higgs-sector observables including higher-order corrections we use the 
soft SUSY-breaking parameters of the best-fit points in the CMSSM and the NUHMl, 
respectively. We show in Fig. [2lthe 5-0 discovery contours for the three decay channels 
H,A — > T+T^ ^ jets (solid lines), jet + ji (dashed lines) and jet + e (dotted lines). The 
parameter regions above and to the left of the curves are within reach of the LHC with 
about 30 fb^ of integrated luminosity. We see that most of the highest-CL regions 
lie beyond this reach, particularly in the CMSSM. At the ILC(IOOO) masses up to 
Ma < 500 GeV can be probed. Within the CMSSM this includes the best-fit point, and 
within the NUHMl nearly the whole 68% C.L. area can be covered. 
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FIGURE 2. The correlations between Ma and tan j3 in the CMSSM ( left panel ) and in the NUHMl ( right 
panel) H^. Also shown are the 5-G discovery contours for observing the heavy MSSM Higgs bosons H,A 
in the three decay channels H,A T^T~ — *■ jets (solid line), jet + /i (dashed line), jet + e (dotted line) 
at the LHC. The discovery contours have been obtained using an analysis that assumed 30 or 60 fb^^ 
collected with the CMS detector ^Z/l [t^/ . 

Finally we discuss the likelihood functions for Mh within the CMSSM and NUHMl 
frameworks obtained when dropping the contribution to from the direct Higgs 
searches at LEP. The results are shown in the left and right panels of Fig. [3l respec- 
tively. The left plot updates that for the CMSSM given in [4]. 

It is well known that the central value of the Higgs mass in a SM fit to the precision 
electroweak data lies below 100 GeV |78, [t^I, but the theoretical (blue band) and 
experimental uncertainties in the SM fit are such that they are still compatible at the 
95% C.L. with the direct lower limit of 1 14.4 GeV IstIi derived from searches at LEP. In 
the case of the CMSSM and NUHMl, one may predict A//, on the basis of the underlying 
model parameters, with a 1-C7 uncertainty of 1 .5 GeV [40], shown as a red band in Fig.[3l 
Also shown in Fig. [3] are the LEP exclusion on a SM Higgs (yellow shading) and the 
ranges that are theoretically inaccessible in the supersymmetric models studied (beige 
shading). The LEP exclusion is directly applicable to the CMSSM, since the h couplings 
are essentially indistinguishable from those of the SM Higgs boson [1381 1391]. but this is 
not necessarily the case in the NUHMl, as discussed earlier. 

In the case of the CMSSM, we see in the left panel of Fig. [3] that the minimum 
of the function occurs below the LEP exclusion limit. While the tension between 
the function for arising from the CMSSM fit and the LEP exclusion limit has 
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FIGURE 3. The functions for Mh in the CMSSM (left) and the NUHMl (right) f^, including the 
theoretical uncertainties (red bands). Also shown is the mass range excluded for a SM-like Higgs boson 
(yellow shading), and the ranges theoretically inaccessible in the supersymmetric models studied. 



slightly increased compared to the earlier analysis performed in [4], the fit result is still 
cornpatible at the 95% C.L. with the search limit, similarly to the SM case. As we found 
in [16|] a global fit including the LEP constraint has acceptable x^. In the case of the 
NUHMl, shown in the right panel of Fig.[3l we see that the minimum of the function 
occurs above the LEP lower limit on the mass of a SM Higgs. Thus, within the NUHMl 
the combination of all other experimental constraints naturally evades the LEP Higgs 
constraints, and no tension between M/, and the experimental bounds exists. 



ACKNOWLEDGMENTS 

We thank O. Buchmiiller, R. Cavanaugh, A. De Roeck, J. Ellis, H. Flacher, G. Isidori, 
K. Olive, F. Ronga and G. Weiglein with whom the results presented here have been 
obtained. This work has been supported in part by the European Community's Marie- 
Curie Research Training Network under contract MRTN-CT-2006-035505 'Tools and 
Precision Calculations for Physics Discoveries at Colliders' (HEPTOOLS). 



REFERENCES 

1. H. P. Nilles, Phys. Rept. 110 (1984) 1. 

2. H. E. Haber and G. L. Kane, Phys. Rept. 117 (1985) 75. 

3. R. Bai-bieri, S. Ferrara and C. A. Savoy, Phys. Lett. B 119 (1982) 343. 

4. O. Buchmueller et ai, Phys. Lett. B 657 (2007) 87 [ai-Xiv:0707.3447 [hep-ph]]. 

5. O. Buchmueller et ai, JHEP 0809 (2008) 1 17 [arXiv:0808.4128 [hep-ph]]. 

6. O. Buchmueller et ai, to appear in Eur. Phys. J. C, arXiv:0907.5568 [hep-ph]. 

7. W. de Boer and C. Sander, Phys. Lett. B 585 (2004) 276 [arXiv:hep-ph/0307049]. 

8. G. Belanger, F. Boudjema, A. Cottrant, A. Pukhov and A. Semenov, Nucl. Phys. B 706 (2005) 41 1 
[arXiv:hep-ph/0407218]. 



9. J. R. ElUs, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Rev. D 69 (2004) 095004 [arXiv:hep- 
ph/03 10356]. 

10. J. R. Ellis, S. Heinemeyer, K. A. Olive and G. Weiglein, JHEP 0502 (2005) 013 [arXiv:hep- 
ph/0411216]. 

11. J. R. ElUs, D. V. Nanopoulos, K. A. Olive and Y. Santoso, Phys. Lett. B 633 (2006) 583 [arXiv:hep- 
ph/0509331]. 

12. J. R. Ellis, S. Heinemeyer, K. A. Olive and G. Weiglein, JHEP 0605 (2006) 005 [arXiv:hep- 
ph/0602220]. 

13. J. ElUs, S. Heinemeyer, K.A. OUve, A.M. Weber, G. Weiglein, JHEP 08 (2007) 083 [arXiv:0706.0652 
[hep-ph]] 

14. J. Ellis, T. Hahn, S. Heinemeyer, K. A. Olive and G. Weiglein, JHEP 0710 (2007) 092 
[arXiv:0709.0098 [hep-ph]]. 

15. J. R. Elhs, S. Heinemeyer, K. A. Olive and G. Weiglein, Phys. Lett. B 653 (2007) 292 
[arXiv:0706.0977 [hep-ph]]. 

16. S. Heinemeyer, X. Miao, S. Su and G. Weiglein, JHEP 0808 (2008) 087 [arXiv:0805.2359 [hep-ph]]. 

17. P Bechde, K. Desch and P. Wienemann, Comput. Phys. Conmiun. 174 (2006) 47 [arXiv:hep- 
ph/0412012]. 

18. R. Lafaye, T. Plehn, M. Rauch and D. Zerwas, Eur. Phys. J. C 54 (2008) 617 [arXiv:0709.3985 
[hep-ph]]. 

19. E. A. Baltz and P Gondolo, JHEP 0410 (2004) 052 [arXiv:hep-ph/0407039]. 

20. B. C. Allanach and C. G. Lester, Phys. Rev. D 73 (2006) 015013 [arXiv:hep-ph/0507283]. 

21. B. C. Allanach, Phys. Lett. B 635 (2006) 123 [arXiv:hep-ph/0601089]. 

22. B. C. Allanach, C. G. Lester and A. M. Weber, JHEP 0612 (2006) 065 [arXiv:hep-ph/0609295]. 

23. B.C. Allanach and C. G. Lester, Comput. Phys. Commun. 179 (2008) 256 [arXiv:0705.0486 [hep- 
ph]]. 

24. B.C. Allanach, K. Cranmer, C. G. Lester and A. M. Weber, JHEP 0708, 023 (2007) [arXiv:0705.0487 
[hep-ph]]. 

25. B.C. Allanach and D. Hooper, JHEP 0810 (2008) 071 [arXiv:0806.1923 [hep-ph]]. 

26. F. Feroz, B. C. Allanach, M. Hobson, S. S. AbdusSalam, R. Trotta and A. M. Weber, JHEP 0810 
(2008) 064 [arXiv:0807.4512 [hep-ph]]. 

27. R. R. de Austri, R. Trotta and L. Roszkowski, JHEP 0605 (2006) 002 [arXiv:hep-ph/0602028]. 

28. L. Roszkowski, R. R. de Austri and R. Trotta, JHEP 0704 (2007) 084 [arXiv:hep-ph/061 1 173]. 

29. L. Roszkowski, R. Ruiz de Austri and R. Trotta, JHEP 0707 (2007) 075 [arXiv:0705.2012 [hep-ph]]. 

30. L. Roszkowski, R. R. de Austri, J. Silk and R. Trotta, Phys. Lett. B 671 (2009) 10 [arXiv:0707.0622 
[astro-ph]]. 

31. R. Trotta et ai, JHEP 0812 (2008) 024 [arXiv:0809.3792 [hep-ph]]. 

32. S. S. AbdusSalam et aL, arXiv:0906.0957 [hep-ph]. 

33. G. Belanger, F. Boudjema, A. Pukhov and R. K. Singh, arXiv:0906.5048 [hep-ph]. 

34. S. S. AbdusSalam et ai, arXiv:0904.2548 [hep-ph]. 

35. P Bechtle, K. Desch, M. Uhlenbrock and R Wienemann, arXiv:0907.2589 [hep-ph]. 

36. LEP Supersymmetty Working Group, lepsusy . web . cern . ch/lepsusy/ . 

37. R. Barate et al. [ALEPH, DELPHI, L3, OPAL Collaborations and LEP Working Group for Higgs 
boson searches], Phys. Lett. B 565 (2003) 61 [arXiv:hep-ex/0306033]. 

38. J. R. Ellis, S. Heinemeyer, K. A. OUve and G. Weiglein, Phys. Lett. B 515 (2001) 348 [arXiv:hep- 
ph/0105061]. 

39. S. Ambrosanio, A. Dedes, S. Heinemeyer, S. Su and G. Weiglein, Nucl. Phys. B 624 (2002) 3 
[arXiv:hep-ph/0 106255]. 

40. G. Degrassi, S. Heinemeyer, W. HoUik, R Slavich and G. Weiglein, Eur. Phys. J. C 28 (2003) 133 
[arXiv:hep-ph/0212020]. 

41. S. Schael et al. [ALEPH, DELPHI, L3, OPAL Collaborations and LEP Working Group for Higgs 
boson searches], Eur. Phys. J. C 47 (2006) 547 [arXiv:hep-ex/0602042]. 

42. P Bechtle, O. Brein, S. Heinemeyer, G. Weiglein and K. E. WilUams, arXiv:081 1.4169 [hep-ph]; 
arXiv:0909.4664 [hep-ph]; see www . ippp . dur . ac . uk/HiggsBounds . 

43. B. C. AUanach, Comput. Phys. Commun. 143 (2002) 305 [arXiv:hep-ph/0104145]. 

44. S. Heinemeyer, W. HolUk and G. Weiglein, Eur. Phys. J. C 9 (1999) 343 [arXiv:hep-ph/9812472]. 



45. S. Heinemeyer, W. HoUik and G. Weiglein, Comput. Phys. Commun. 124 (2000) 76 [arXiv:hep- 
ph/9812320]. See www . f eynhiggs . de . 

46. M. Frank, T. Hahn, S. Heinemeyer, W. HolUk, H. Rzehak and G. Weiglein, JHEP 0702 (2007) 047 
[arXiv:hep-ph/06 11326]. 

47. T. Moroi, Phys. Rev. D 53 (1996) 6565 [Erratum-ibid. D 56 (1997) 4424] [arXiv:hep-ph/9512396]. 

48. G. Degrassi and G. F. Giudice, Phys. Rev. D 58 (1998) 053007 [arXiv:hep-ph/9803384]. 

49. S. Heinemeyer, D. Stockinger and G. Weiglein, Nucl. Phys. B 690 (2004) 62 [arXiv:hep- 
ph/03 12264]. 

50. S. Heinemeyer, D. Stockinger and G. Weiglein, Nucl. Phys. B 699 (2004) 103 [arXiv:hep- 
ph/0405255]. 

51. M. Davier, Nucl. Phys. Proc. Suppl. 169 (2007) 288 [arXiv:hep-ph/0701 163]. 

52. J. P Miller, E. de Rafael and B. L. Roberts, Rept. Prog. Phys. 70, 795 (2007) [arXiv:hep-ph/0703049]. 

53. F Jegerlehner, Acta Phys. Polon. B 38, 3021 (2007) [arXiv:hep-ph/0703125]. 

54. M. Passera, W. J. Marciano and A. Sirhn, Phys. Rev. D 78, 013009 (2008) [arXiv:0804. 1 142 [hep- 
ph]]. 

55. M. Davier, A. Hoecker, B. Malaescu, C. Z. Yuan and Z. Zhang, arXiv:0908.4300 [hep-ph]. 

56. M. Davier et aL, arXiv:0906.5443 [hep-ph] . 

57. G. Isidori and P Paradisi, Phys. Lett. B 639 (2006) 499 [arXiv:hep-ph/0605012]. 

58. G. Isidori, E Mescia, P Paradisi and D. Temes, Phys. Rev. D 75 (2007) 115019 [arXiv:hep- 
ph/0703035], and references therein. 

59. F. Mahmoudi, Comput. Phys. Conmiun. 178 (2008) 745 [arXiv:0710.2067 [hep-ph]] and 
arXiv:0808.3 144 [hep-ph]. 

60. D. Eriksson, E Mahmoudi and O. Stal, JHEP 0811 (2008) 035 [arXiv:0808.3551 [hep-ph]]. 

61. S. Heinemeyer, W. Holhk, D. Stockinger, A. M. Weber and G. Weiglein, JHEP 0608 (2006) 052 
[arXiv:hep-ph/0604147]. 

62. S. Heinemeyer, W. Holhk, A. M. Weber and G. Weiglein, JHEP 0804 (2008) 039 [arXiv:0710.2972 
[hep-ph]]. 

63. G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 176 (2007) 367 
[arXiv:hep-ph/0607059]. 

64. G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 149 (2002) 103 
[arXiv:hep-ph/0 112278]. 

65. G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 174 (2006) 577 
[arXiv:hep-ph/0405253]. 

66. P Gondolo et ai. New Astron. Rev. 49 (2005) 149. 

67. P Gondolo et aL, JCAP 0407 (2004) 008 [arXiv:astro-ph/0406204J. 

68. P Skands et aL, JHEP 0407 (2004) 036 [arXiv;hep-ph/031 1 123]. 

69. B. Allanach et aL, Comput. Phys. Conomun. 180 (2009) 8 [arXiv:0801.0045 [hep-ph]]. 

70. G. Aad et aL [The ATLAS Collaboration], arXiv:0901.0512. 

71. G. L. Bayatian et aL, CMS Collaboration, J. Phys. G34, 995 (2007); see: 
cms doc . cern . ch/ cms/ cpt/tdr/ . 

72. K Gianotti et aL, Eur. Phys. J. C 39 (2005) 293 [arXiv:hep-ph/0204087]. 

73. TESLA Technical Design Report [TESLA CoUaboration], arXiv:hep-ph/01063 15; see: 
tesla . desy . de/ new_pages/TDR_CD/ start .html . 

74. J. Brau et aL [ILC Collaboration], ILC Reference Design Report Volume 1 - Executive Summary, 
arXiv:0712.1950 [physics.acc-ph]; 

G. Aarons etaL [ILC Collaboration], International Linear Collider Reference Design Report Volume 
2: Physics at the ILC, arXiv:0709.1893 [hep-ph]. 

75. M. S. Carena, S. Heinemeyer, C. E. M. Wagner and G. Weiglein, Eur. Phys. J. C 26 (2003) 601 
[arXiv:hep-ph/0202167|. 

76. M. S. Carena, S. Heinemeyer, C. E. M. Wagner and G. Weiglein, Eur. Phys. J. C 45 (2006) 797 
[arXiv:hep-ph/0511023]. 

77. S. Gennai, S. Heinemeyer, A. Kalinowski, R. Kinnunen, S. Lehti, A. Nikitenko and G. Weiglein, Eur. 
Phys. J. C 52 (2007) 383 [arXiv:0704.0619 [hep-ph]]. 

78. M. Goebel, talk given at Rencontres de Moriond EW2009 . 

79. LEP Electroweak Working Group, lepewwg . web . cern . ch/LEPEWWG/Welcome . html . 



